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Abstract 

The physiological and morphological traits of flag leaf play an important role in determining the biomass and seed yield of 
crop. In order to understand the genetic control of different traits of flag leaf, a recombinant inbred lines (RIL) population 
derived from a cross between ‘Roshan’ and ‘Falat’ (seri82) was used to identify the quantitative trait locus (QTL) underlying 
4 morphological traits under control and saline environments. We used QTL Cartographer and QTL Network to estimate QTLs 
of phenotypic values and 16 additive QTLs for different traits of flag leaf were detected by using a genetic map of 730 DArT 
and SSR markers. We also identified one pairs of epistatic effects using the phenotypic value of single environment for one 
trait (flag leaf width) in 2 chromosomes. The QTL Network is used to detect QTLs and interaction between them and | additive 
QTL for flag leaf length and 1 additive QTL for flag leaf area were identified. None of the additive QTLs had interaction with 
environment. Studies of the locus linked to the flag leaf length and weight of flag leaf QTLs revealed that wPt-0549 is belonged 
to the part of sequence that could encode serpin and cytochrome P450 protein. It could be concluded that serpin and cytochrome 
P450 genes which are consist of wPt-0549 regulate leaf senescing through their regulation under salinity stress. 

Keywords: Flag leaf traits; Wheat; Quantitative trait loci; Salt tolerance 

Abbreviations: RIL: recombinant inbred lines, QTL: quantitative trait locus, EC: electrical conductivity, FLL and FLW: flag 
leaf length and width of the flag leaf, FLA: flag leaves area, WFL: weight of flag leaves, ANOVA; Analysis of variance, SPSS: 








Statistics 15 software, SSRs: simple sequenced repeats, CIM: composite interval mapping, CV: coefficient of variation 
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Introduction 


Improvements in yield factors (such as biological yield, 
thousand grain weight, kernel per spike, spike weight) and its 
related traits can make seed yield better [1]. The inheritance 
of yield is not simple and not only associated to kernel weight, 
kernel number per spike and spike number per unit area but 
also influence by additional things, for example the storage 
and transportation of photosynthetic yields [2]. It is appeared 
that plant parts under the flag leaf node have a small role in 
final performance of plant. The photosynthetic potential 
which impressing the yield and stress reaction of wheat plants 
was affected by leaf size, shape, posture as well as ear and 
awn [3, 4]. As a result, it seems logical that plant structures 
above the flag leaf node define the final performance of the 
plant [5]. Seed yield affect by the leaves, especially the flag 
leaf, where the most of photosynthetic activities (almost 50 
%) controlled by it [6, 7]. Grain filling in cereals is usually 
done by three upper leaves, which are the primary source of 
photosynthesis in plant [8, 9, 10]. Flag leaf in wheat has an 
important role in the proportion of photosynthetic assimilates 
compared to other leaves and can produce a large portion of 
the carbohydrates stored in the seed [10-12]. So flag leaf in 
wheat is considered the most important component in seed 
yield [13-14]. The importance of flag leaf is because of the 
short interval between flag leaf and spike and the other is that 
flag leaf stays green longer than other leaves [15]. About 30 
to 50 percent of dry matter at maturity is related to flag leaf, 
and it is the most important site for photosynthesis in the seed 
filling stage [16, 17]. 

The performance of a plant is a relationship between many of 
the physiological and biochemical processes occurring in 
plants, and a number of genes are involved in this trait and the 
environment can affect them [18, 19]. Since flag leaf plays an 
important role, its size is probably important [6]. 

The positive effect of leaf size on 1000-grain weight, number 
of grains per spike, grain weight per plant or other yield- 
related traits has been confirmed by several studies in cereals 
like rice and wheat [2, 20-22]. By improving the flag leaf 
morphology, photosynthetic efficiency and yield can be 


expected to improve [23]. 


Many researchers observed that decrease in leaf area caused 





a reduction in plant growth and beside that plant development 
and leaf area decrease by salinity [24]. Salt stress reduce leaf 
area and plant height [25]. Death of leaves before 
development can result from a high level of salinity [26]. The 
photosynthetic efficiency of next transpiring leaves like flag 
leaf was affected by this event in adult wheat plants [26]. 

A new vision for crop growth and development can be 
achieved by understanding the details function of 
physiological and morphological characters of flag leaf on 
plants performance. The analysis of complex traits like yield 
can be performed using the genetic maps and molecular 
markers for continuous distribution expression of 
phenotypes. They are also applicable in detecting related loci 
and measuring the average genes involved in a trait and 
defining the relations among QTLs and the environment by 
using quantitative genetics relied upon statistics [27]. One of 
the factors that affects quantitative traits is the relation 
between QTLs and the environment [28-31]. A number of 
loci with an effect in flag leaf traits have been identified in 
rice, barley and wheat [19, 20, 23, 32-39]. However, more 
research is needed to identify QTLs for flag leaf morphology, 
especially in different environments. 

In current study for detecting loci affecting morphological 
traits of flag leaf at flowering stage by using approach of QTL 
Cartographer and QTLNetwork [40], we studied a 
recombinant inbred lines (RIL) population derived from a 
cross between ‘Roshan’ and ‘Falat’ (seri82) in normal and 
saline environments. The aim and significance of this study 
in both theory and application was to study the relationships 
among flag leaf traits, the effect of salinity on flag leaf traits 
and to use identified QTLs in marker assisted selection for 
improving plant yield by improving flag leaf traits, for 


example by increasing the photosynthetic area. 


Materials and Methods 


Plant materials and field trials 

In present study we use 319 inbred lines (F7) derived from a 
cross between Roshan and Falat cultivar. The Roshan is a salt 
tolerant cultivar and the Falat (seri82) is a salt sensitive 
cultivar [41]. These two parents have an apparently 


significant difference in terms of morphological traits related 
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to flag leaf. 

The population with their parents were grown in two 
environments; Yazd (55 31°N 22 54’E, 1237 m above sea 
level, as control environment) and Ardakan (53 48°N 32 
20°E, 1237 m above sea level, as saline environment), both 
located in the central part of Iran. The soil EC of Yazd was 
4.5 dS m'! and the water EC was 3.5 dS m’!. Soil samples 
were collected from five representative locations (upper left, 
upper right, lower left, lower right and the middle) of the 
fields in both Yazd and Ardakan for electrical conductivity 
(EC) measuring [42]. The soil EC of Ardakan was 12 dS m'! 
and the water EC was 10 dS m'!. The experiment was 
performed in an alpha lattice design with two replicates in 
each environment. The material was sown in two 2 m rows 
(interplant distance 2 cm, inter-row distance 20 cm). At 
flowering, flag leaf length and width of the flag leaf (FLL and 
FLW) were measured using 10 samples from each plot in both 
environments. The flag leaf length and width were calculated 
as follows; flag leaf length: measuring the leaves from the 
beginning of the ligula to the end of the tip of the leaf; and 
flag leaf width: measuring the leaves from the widest part of 
the flag leaf. To find out the flag leaves area (FLA), leaf 
length and width based on centimeter was measured and 
multiplied by the coefficient 0.74 [43]. We oven-dried the 
flag leaves at 75 ‘C for 48h and measured the weight of flag 
leaves (WFL) in both environment. Agronomic and 
phenological traits like seed yield, biological yield, thousand 
seed weight, plant height, spike length, seed weight per spike, 
number of seed per spike and harvest index were also 
recorded at maturity using 20 random plants in each plot. In 
order to measure seed yield, each plot was completely 


harvested. 


Analysis of variance (ANOVA) of each data set was done 
using SPSS Statistics 15 software (SPSS Inc, Chicago, USA). 





Data and QTL analysis 


The average of two replicates for each trait was calculated 
before QTL analysis. For QTL analysis, the linkage map of 
the ‘‘Roshan = Falat’’ population was used [44]. This map 
contained all the wheat chromosomes and includes 730 
markers, including 705 DArT and 25 simple sequenced 
repeats (SSRs). We determined the main effect QTL for each 
single-environment phenotypic value by composite interval 
mapping (CIM) analysis using Win QTL Cartographer and 
QTLNetwork software and also pleiotropic QTLs by multi- 
traits composite interval mapping (MCIM) method [45] using 
Win QTL Cartographer. We also determined the effect of a 
and aa epistatic QTL and their environmental interactions (ae 
and aae) with QTL Network software version 2.1 [46, 47] 
with a mixed linear composite interval mapping method and 
with joint analysis of the multi-environment phenotypic 
values for FLL, FLW, FLA and WEL. In both software, 
significant LOD thresholds were determined by 1,000 
permutation test and p = 0.05. 

After QTL mapping, in order to reveal the relationship 
between the identified QTL for leaf traits with a strong LOD 
and the role of the genes, the flag leaf traits locus and 
characteristics of their related genes were studied using 


ensemble plant (link) and NCBI [48]. 


Results 

Evaluation of the flag leaf phenotype and correlations 
among the traits 

The salt-tolerant Roshan cultivar had higher phenotypic 
values than sensitive Falat in both environments (Fig 1). 
There were significant differences in all traits between both 
parents and inbred lines in both environments (Supplemental 


Table 1). 
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Figure 1: Frequency distribution of the means of traits and parents positions among 319 recombinant inbred lines (RIL) under stress environments. 
R and F are Roshan and Falat, respectively. 
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Supplemental Table 1: Phenotypic performance of flag leaf traits in two parents and RIL population under Yazd (control) 


and Ardakan (saline) environments. 



































Trait environment parent RILs 
Roshan Falat Mean Min. Max. SD CV (%) 
* 
Flag leaf length, Yazd 20 + 1.27 16.03 + 2.53 15.19 9.43 22.2 2.29 15.09 
FLL (cm) Ardakan 15.5 + 0.26 10.8 + 2* 10.89 6.88 16.03 1.48 13.62 
* 
Flag leaf width, Yazd 14.741.2 12.841 12.73 8.61 17.03 1.57 12.33 
FLW (mm) Ardakan 14.541 12+1* 11.91 8.17 16.17 1.38 11.6 
Yazd 1895.66 + 122.23 1283.16 + 60.55* 1161.82 633.58 1996.6 239.71 20.63 
Flag leaf area, FLA Ardakan 1457.51 +25.41 845.97 + 208.8* 832.13 410.81 1300.21 159.17 19.13 
Yazd 0.185 + 0.005 0.144 + 0.008* 0.12 0.08 0.19 0.03 23.2 
Weight of flag leaf, = 
WFL Ardakan 0.138 + 0.004 0.099 + 0.007* 0.08 0.04 0.14 0.02 25.) 





























“The parents were significantly different at the 0.05 probability level 


The phenotypic values for the flag leaf traits exhibited wide 
ranges between the 319 RILs, and the coefficient of variation 
(CV) of all traits were of 15 % except for FLL and FLW in 
normal environment (Yazd environment: the N environment) 


FLL and FLA _ displayed 


(Supplemental 


Table 1). 


environment effects. The mean values of FLL, FLW and FLA 
among RILs were higher in N environment than that in the S 
environment (Ardakan environment: the S environment); 
however, this difference for FLW was not significant. We 


observed significant transgressive segregation in both 





www.acquirepublications.org/JARL 





Journal of Agricultural Research and Livestock 


directions and continuous variation for the frequency 
distribution among all measured traits due to their polygenic 
inheritance patterns (Fig 1). The significant phenotypic 
correlations were observed among flag leaf traits and several 
measured traits in normal and stress environments 
(Supplemental Table 2). We detected a positive and 
significant correlation between FLL and FLW in both 
environments. 


The result showed that there were a positive and significant 


correlation between FLL, FLW and FLA with spike length, 





number of fertile spikelets per spike and total spikelet per 
spike in both environments 

Also, WFL had a positive and significant correlation with 
most of the yield components. 

A positive correlation was observed between seed yield and 
leaf area in stress environment. Also, a significant correlation 
was observed between FLW and seed yield in stress 


environment. 


Supplemental Table 2: Correlation coefficients between flag leaf traits and other traits under Yazd (above) and Ardakan 


(below) environments. 













































































SY! | BY? | TSW? | PH‘ SL‘ SW® | SW\S’ | NS\S® | AL? | SS\S!° | TS\S! | FS\S? | HI? | FLA™ | FLW! 
FLL” | -0.08 | 0.11° | 0.01 0.06 | 0.25" | 0.02 0 -0.01 | 0.28" | 0.14" | 0.26" | 0.19" | -0.14" | 0.83" | 0.19% 
FLW -0.06 | 0.02 0.01 0.03 | 0.13" | 0.18 | 0.17 | 0.15" | -0.11" | 0.15™ | 0.31" | 0.24 | -0.08 | 0.70% 
FLA -0.09 | 0.09 0.01 0.06 | 0.25" | 0.11* 0.09 0.08 | 0.14 | 0.19 | 0.36" | 0.27™ 0 15* 0.68" 
WFL 0.04 | 0.15" | -0.03 | 0.08 | 0.31 | 0.15" ) 0.13" | 0.18"* | 0.11° | 0.16 | 0.30°* | 0.25™ 0 19°" 0.76" | 0.79% 
FLL 0.04 | 0.07 | 0.14" | 0.08 | 0.24" | O01 0.08 0.07 | 0.20" | 0.01 0.26" | 0.25" | -0.04 | 0.80 | 0.15™ 
FLW 0.1 | 0.18" | -0.03 | 0.11" | 0.23" | 0.22™ | 0.20" | 0.18" | -0.1 -0.02 | 0.27 | 0.28" | -0.07 | 0.70" 
FLA 011° | 0.16" | 0.08 | 0.13" | 0.31 | 0.21™ | 0.18" | 0.17" | 0.09 | -0.01 | 0.35 | 0.35" | -0.07 0.70" 
WFL 0.02 | 0.13" | -0.06 0.1 | 0.27 | 0.22™ | 0.19" | 0.23" | 0.02 0.02 0.31" | 0.30" 0 1s” 0.73" | 0.76" 





“Indicates significance at the 0.05 probability level 


“Indicates significance at the 0.01 probability level 


' Seed yield; 7 biological yield; + thousand seed weight; + plant height; > spike length; ’ seed weight per spike; ® number of seed per spike; ? awn 


length; !° sterile spikelet per spike; !! total spikelet per spike; !* fertile spikelet per spike; !* harvest index; '* flag leaf area; ! flag leaf width; !° 


flag leaf K*+/Na* concentration ratio; !7 flag leaf K+ concentration; '* flag leaf Na* concentration; ! flag leaf greenness; 7° weight of flag leaf; 7! 


flag leaf symptoms; ”” flag leaf length 


Mapping QTLs with additive effects for flag leaf traits 
QTL analysis of the 313 RILs detected 16 additive QTL on 9 
chromosomes for 4 traits by using Win QTL Cartographer 
and also 3 additive QTLs on 3 chromosomes for 2 traits using 
QTLNetwork software (Table 1 and 2; Fig 2). 

The distribution of identified QTLs in A, B and D genomes 
was 4,13 and 2 respectively. The 6 of the QTLs identified by 


QTL Cartographer software were found in the Yazd 
environment while the remaining 10 were manifest in the 
Ardakan environment. The QTLs identified by QTLNetwork 
software were found only in Ardakan environment. The QTL 
explained 3.77 to 12.31 % of the phenotypic variation (Table 
1 and 2). 
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Table 1: Location of additive QTLs flag leaf traits in RIL population under Yazd (control) and Ardakan (saline) 


environments by CIM 


























h(a 
Traits ae ees site» | Range a Marker | citee | Range | environments | a or aa or aael aae. | h?(aae) 
i interval (cM) j interval (cM) a) 
wPt- 
FLA- | 743486- 5.6- F : 2 : _ « ‘ 
FLA TBI wPt- 10.4 147 Ardakan 46.3705 8.4 
8283 
tPt- 
WFL- 3689- 53.1- 4 7 : . 2 we x ” 
WFL 6B1 wPt- 58.8 68.4 Ardakan 0.0046 3.9 
8554 
wPt- 
WFL- | 742051- 164.1- |, : : ; * r : 
7A] Pe 171.1 178.4 Ardakan 0.0081 8.6 
3226 
wPt- wPt- 
FLW- 5906- 120.5- | FLW- 1761- 41.1- saa * * 
FLW 3B1 wPt- 126.1 B21 6B1 wPt- 49.1 55.1 Ardakan -0.504 8.2 
1171 4924 
wPt- 
FLL on 5730- | 15.7 ae * * * * Both 4.3268! 107 5 ge oan | 0.071 
tPt-5597 : : 
wPt- 
FLA - | 743486- 5.6- « : - ai 7 
FLA BI wPt- 10.4 147 Both 47.9176 3.32 0 0 0.0001 
8283 















































> Site means the distance of F value peak for QTL after the first marker in the marker interval 


Table 2: QTLs with additive effects (a) (above the highlight parts), epistatic effects (aa) (the highlight parts) and with additive 


effects (a) and additive x environment interaction effects (ae) (below the highlight parts) under Yazd (e7) and Ardakan (ez) 


environments detected by QTLNetwork. 













































































Traits Chr.? Marker interval Range (cM) Environments LOD R? (%) Add 
FLL 2B-1 wPt-3561-wPt-8072 12.2-17.2 Yazd 3.7351 4.71 0.5193 
7B-1 wPt-2273-wPt-4814 64.7-92.1 Yazd 3.2948 5.57 -0.5517 
3A-1 wPt-6422-wPt-0549 37.2-68.9 Ardakan 3.2912 12.31 0.5337 
3B-2 wPt-4412-wPt-7614 14.8-20.7 Ardakan 4.3665 5.06 -0.3873 
5D rPt-3825-wPt-667413 82.3-96.4 Ardakan 3.9951 5.28 0.4083 
FLW 2B-1 wPt-741382-wPt-0047 105.4-135.8 Ardakan 2.8708 4.71 0.3219 
3B-1 wPt-11295-Xgwm566 388-406.1 Ardakan 3.1279 5.58 -0.343 
6B-1 wPt-74323 1-wPt-2498 47.9-68.5 Ardakan 3.1859 4.56 0.339 
FLA 2B-1 wPt-3561-wPt-8072 11.9-17.2 Yazd 3.6184 4.31 53.56 
4D wPt-666459-wPt-672143 2.3-10.7 Yazd 3.6429 5.22 -64.9438 
6B-1 tPt-3689-wPt-8554 55.8-65.4 Yazd 3.3475 4.1 52.4582 
7B-1 wPt-743486-wPt-8283 6.3-16.5 Yazd 3.1884 3.77 -47.5711 
TA-1 wPt-3226-wPt-2525 178-181.4 Ardakan 3.3082 3.88 -35.2663 
7B-1 wPt-743486-wPt-8283 7.1-16 Ardakan 4.2814 5.06 -36.8172 
WFL 3A-1 wPt-0549-wPt-9154 50.2-81.4 Ardakan 2.7151 11.82 0.0072 
6B-1 wPt-74323 1-wPt-2498 37.5-52 Ardakan 3.3394 4.43 0.0046 





ns 
2 


Not significant and significant at 0.001 probability level, respectively 


@ Marker interval means the interval of the F value peak for QTLs 
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A total of 3 QTL on 3 chromosomes were identified for flag 
leaf width which were contributed by parents’ alleles. The B 
genome had the largest number of QTLs (All QTL) for this 
trait. All of the QTLs were identified by QTL Cartographer 
software. 

Of these QTLs for FLW, one which was located on 
chromosome 6B - 1(wPt-743231-wPt-2498), had the same 
location with a QTL affecting WFL in the Ardakan 
environment and explained 4.56 % of the phenotypic 
variation. Both QTLs were contributed by ‘Roshan’ alleles. 
Of the 3 putative QTLs for this trait, all of them were 
identified under salt stress environments (Ardakan). Both 
parental alleles contributed to the effect of flag leaf width 
under the salt stress environment. 

This study identified seven QTL on 5 chromosomes for the 
flag leaf area which were contributed by parents’ alleles. The 
B genome had the largest number of QTLs (5 QTL) for this 
trait. From these 7 QTLs, four of them were identified by 
QTL Cartographer software and only one on chromosome 
7B-1 was identified in both environments and also it was 
identified in salt stress environment by QTLNetwork 
software. 

One of these QTLs located on chromosome 2B - | (wPt-3561- 
wPt-8072) had the same location with a QTL for FLL in the 
Yazd environment and explained 4.31 % of the phenotypic 
variation. Both contributed by ‘Roshan’ alleles. Of the 6 
putative QTLs identified by QTL Cartographer for this trait, 
one of them on chromosome 7B - | (wPt-743486-wPt-8283) 
was identified at the same position under normal and salt 
stress environments. Under normal environment, two QTLs 
for flag leaf area from Falat had a negative effect. For the 
others in this environment, Roshan alleles were associated 
with positive effects on flag leaf area. The Falat alleles 
contributed to the effect of flag leaf area under salt stress 
environment. 

This study detected four QTL on 3 chromosomes for the flag 
leaf weight which were contributed by parents’ alleles. Both 
the A and B genome had the same number of QTLs (2 QTL) 
for this trait. From these 4 QTL, two of them were identified 
by QTL Cartographer software and the rest were identified by 
QTL Network software. 


One of them on chromosome 3A - 1 (wPt-0549-wPt-9154) 
was co-located with a QTL affecting FLL in the Ardakan 





environment and explained 11.82 % of the phenotypic 
variation and the other one on chromosome 6B - | (wPt- 
74323 1-wPt-2498) was co-located with a QTL affecting 
FLW. Both contributed by ‘Roshan’ alleles. 

Blast result indicated that the sequence of a DArt marker 
(wPt-0549) which was linked to flag leaf length and weight 
of flag leaf QTLs in Ardakan environment in chromosome 3A 
- 1 showed this sequence located in 5B chromosome of 
Triticum aestivum and overlap with TraesCS5B02G48 1400 
gene with 99.5 % _ of identity. According to Ensembl Plants 
results, wPt-0549 region could encode coding protein 
belonged to serpin family. Furthermore, we have also 
investigated on NCBI to find any similar sequences using 
blastn. Results showed that wPt-0549 is similar to 
HE774676.1 (contig ctg447) locus of Triticum aestivum 
located at the short arm of 3D chromosome with 98.98 % of 
identity. HE774676.1 locus could encode cytochrome P450 
protein. 

Mapping QTLs with different effects 

To breed a physiologically efficient and productive wheat 
genotype, it would be helpful to know the epistatic gene 
actions operating in the inheritance of the physiological traits, 
such as flag leaf area [1]. We detected a pairs of epistatic 
effects using single-environment phenotypic value for one 
trait on 2 chromosomes (Table 2; Fig. 2) by composite 
interval mapping (CIM) analysis using QTLNetwork 
software. The single epistasis was identified for FLW. The 
FLW-3B1/FLW-6B1 epistasis resulted in decrease FLW in 
recombination type qlq1Q2Q2. 

We also identified 2 QTLs for two traits on two chromosomes 
by the combined analysis of the multi-environment 
phenotypic values under normal and saline conditions, Both 
of the QTLs did not show interacting effects in the 
environment, explaining the phenotypic variation ranging 
between 1.07 to 3.32 % (Table 2). We also identified no pairs 
of epistatic effects for these traits. 

We found one QTL for FLL on chromosome 4A - 1, which 
‘Falat’ alleles contributed to this QTL. Also, we found only 
one QTL for FLA on chromosome 7B - | which ‘Roshan’ 


alleles contributed to this QTL 
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Figure 2: Location of QTLs for flag leaf traits on linkage groups. 
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Detecting QTLs cluster and pleiotropy 

MCIM analysis using the individual environment data 
indicated 4 intervals on chromosomes 2B - 1, 3A - 1, 6B - 1 
and 7B - 1 for additive QTLs. The interval wPt-356/-wPt- 
8072 on chromosome 2B - 1 affects FLL (N) and FLA (N), 
the interval wPt-0549-wPt-9154 on chromosome 3A - 1 
affects FLL (S) and WEL (S), the interval wPt-743231-wPt- 
2498 on chromosome 6B - | affects FLW (S) and WFL (S), 
and the interval wPt-743486-wPt-8283 on chromosome 7B - 
1 affects FLA (N) and FLA (S) (one trait in the both 


treatments). 


Discussion 

Relationships among the flag leaf traits and salt tolerance 
Different physiological mechanisms play a role in salt 
tolerance and are controlled at different stages of the growth 
cycle [49-51]. A complete list of traits and methods used to 
evaluate salt tolerance has been reviewed elsewhere [52]. 
Flag leaf has an important role in cereals like wheat, because 
it provides the maximum amount of photosynthesis 
assimilates to be stored in the grains. It seems that flag leaf 
with higher length had higher width. Also, we found a 
positive and significant correlation between FLA with FLL 
and FLW in both environments, which was in good 
agreement with [53]. Other researchers found similar 
relationship between FLL and FLW [6, 38, 53] . A greater 
flag leaf area will eventually help to increase the 
photosynthetic efficiency by increasing the production of 
photosynthesis, which is then translocated into grains and 
increase their weight. Therefore, flag leaf area has a direct 
relationship to seed yield [54]. As indicated in the results, 
there was a significant correlation between FLW and seed 
yield in stress environment, which was in good agreement 
with [6]. 

[25] reported that salinity reduces the leaf area. The 
morphological traits of flag leaves, including length, width 
and area are inherited quantitatively and influenced greatly by 
the environments [55, 33-34] A positive and significant 
correlation was found between FLA and plant height in stress 
environment, which was in good agreement with other studies 
[15, 56, 57]. Furthermore, other researchers found similar 


relationships between FLL, FLW and FLA with spike length, 


number of fertile spikelets per spike and total spikelet per 





spike similar to our results [15, 57]. 

Roshan is a salt tolerant cultivar and its flag leaf traits include 
FLL, FLW, FLA and WFL were higher than Falat. 

However, in this study there was a positive relationship 
between seed yield and leaf area in the stress environment. 
This positive and significant association of flag leaf area with 
seed yield suggested the importance of flag leaf area in 
evolving superior genotypes [57]. It has been proved that the 
flag leaf, stem and head are the closest source to the grain. 
The flag leaf could produce a large proportion of the 
carbohydrates stored in grains [10-12]. Therefore, flag leaf is 
one of the important components in determining seed yield 
potential in cereal and plays a major role in increasing 
productivity [14]. Physiological studies of wheat have shown 
that flag leaf contribution towards grain weight accounts for 
41- 43 % of dry matter in the kernel at maturity and is the 
major photosynthetic site during the grain filling stage [58]. 
Several researchers observed a positive correlation between 
seed yield and FLA [15, 54, 53, 59, 60]. Therefore, the flag 
leaf area has a direct relationship to seed yield [54]. 
Relationships between flag leaf QTLs 

The analysis of the QTL controlling morphological traits of 
flag leaf has been intensively studied before in some cereal 
[19, 20, 23, 33, 34, 37, 61-65]. However, genetic control of 
flag leaf traits in saline environment has been reported 
somewhat in wheat up to date. 

[36, 23, 38, 66] observed a QTL controlling FLW (flag leaf 
width) on chromosome 5A, [12] observed 3 QTLs controlling 
FLW on chromosomes 1B and 4B, [67] detected eight QTLs 
for FLW on 1A, 1B, 3A, 4D, 5A, 5D, and 7D, [11] detected 
4 QTLs on chromosomes 1A, 1B, 5A and 7A for FLW,. Also 
In barley, there is a major FLW QTL in chromosome 5H 
[34,68]. In the present study, we detected 3 QTLs on 
chromosomes 2B - 1, 3B - | and 6B - | for flag leaf width. 
The roles and relations between QTLs are influenced by 
epistatic effects and pleiotropy; As a result, one minor QTL 
may have a great impact on regulatory pathways [69]. 

[38] detected 7 chromosome regions for flag leaf length on 
chromosomes 1B, 2B, 2D,3B, 4A, 5B, 6D in wheat, [7] 
detected three QTLs on chromosomes2D, 4D, and 5B for 
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FLL, [67] detected three QTLs for FLL on 1A, 5B, and 6A, 
[68] detected 3 QTLs for flag leaf length on chromosomes 

2 QTLs for flag leaf length on the chromosomes 5H and 7H. 
In the present study, we also detected 5 QTLs on 
chromosomes 2B - 1, 3A - 1, 3B - 2, 5D and 7B - 1 for flag 
leaf length. 

There is a strong need to understand the extent of the epistasis 
and QTL x environment (QE) interactions in QTL studies. 
Several studies have shown the interaction between two or 
more QTLs and/ or between QTL and the environment as well 
as the identification of their additive effects [71 - 72]. 
However, there are few reports on the identification of the 
epistatic and Q x E interactions under salt stress conditions, 
especially for the flag leaf traits in wheat. To date, wheat 
researchers have identified different QTLs with epistatic and 
Q x E effects for both physiological and biomass traits [71]. 
In the present study we could find three pairs of epistatic 
effects using single-environment phenotypic value for three 
traits. 

QTL co-location and traits correlation 

QTL clusters or co-located QTLs for different traits in salinity 
conditions have been reported in previous studies [50,72]. 
Masoudi et al. (2015) have reported six cluster on 
chromosomes 1B, 2B (two), 3B (two) and 6B which are at 
least relevant with two different traits. Also co-located QTL 
for flag leaf-related traits have been observed by previous 
researchers in wheat [19,38,39]. 

In the present study, three intervals were identified on 
chromosomes 2B - 1, 3A - | and 6B - 1, which had an effect 
on more than one trait. For each of the three loci, the 
relationships between related traits in these clusters were 
consistent with the additive effects of the corresponding 
QTLs. All three intervals had a similar additive effect for the 
two relevant traits, and the correlation coefficients between 
most of these traits were significant and positive. As already 
mentioned, each cluster may be a single locus or several loci 
linked tightly together [50,71,74]. Four chromosome 
intervals which either affected more than one _ trait 
(pleiotropy) or affected a trait in both N and S environments 
were identified by MCIM analysis. These co-localized QTLs 
were located on chromosomes 2B - 1, 3A - 1, 6B - 1 and 7B 


- 1, and they may be single loci that affect more than one trait. 


Molecular aspects of locus link to flag leaf traits QTLs 





We have observed that wPt-0549 is belonged to the part of 
sequence that could encode serpin protein. It has been shown 
that serpin family is responsive to environmental stress. 
Serpin gene showed up-regulation under salt stress in root 
tissues [75]. On the other hand, It has been also shown that 
Serpin gene interact with RD21 which induced by salt stress 
in senescing leaves and stressed leaf tissues [76-80]. 
Furthermore, we have also observed that wPt-0549 is also 
related to cytochrome P450 protein. It has been shown that 
cytochrome P450 genes were up-regulated during wheat flag 
leaf senescence [81]. Our results indicated that wPt-0549 
marker which is a part of serpin and cytochrome P450 genes 
region showed a significant correlation with flag leaf 
characteristics under salinity stress condition. It has been 
shown that these two genes are involved in leaf senescing. 
Therefore, it could be concluded that serpin and cytochrome 
P450 genes (which are consist of wPt-0549) induce leaf 


senescing through their up-regulation under salinity stress. 


Conclusions 

In this study, using a high-density DArT and SSR linkage 
map, 16 additive QTLs were detected in Yazd and Ardakan 
to control the morphology of wheat flag leaf. One stable QTL 
identified for flag leaf area in the interval wPt-743486-wPt- 
8283 on chromosome 7B - 1. The results showed that most of 
the detected QTLs in the two conditions did not have 
interactions with the environment. In this study, a pairs of 
QTLs with epistatic effects (AA) were identified for flag leaf 
width. 
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